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Abstract 

An analytic solution of the Regge- Wheeler (RW) equation has been found via the 
Frobenius method at the regular singularity of the horizon 2M, in the form of a time 
and radial coordinate dependent series. The RW partial differential equation, de- 
rived from the Einstein field equations, represents the first order perturbations of the 
Schwarzschild metric. The known solutions are numerical in time domain or approxi- 
mate and asymptotic for low or high frequencies in Fourier domain. The former is of 
scarce relevance for comprehension of the geodesic equations for a body in the black 
hole field, while the latter is mainly useful for the description of the emitted gravita- 
tional radiation. Instead a time domain solution is essential for the determination of 
radiation reaction of the falling particle into the black hole, i.e. the influence of the 
emitted radiation on the motion of the perturbing mass in the black hole field. To 
this end, a semi-analytic solution of the inhomogeneous RW equation with the source 
term (Regge- Wheeler-Zerilli equation) shall be the next development. 
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1 Introduction 

Regge and Wheleer M proved the vacuum stability of the Schwarzschild black 
hole, while Zerilli || studied the emitted radiation via the radial wave equa- 
tion for polar perturbations, which source is a freely falling test mass m towards 
the black hole of large mass M. Moncrief B showed the gauge invariant signif- 
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icance of the wave equations. 



Past work was concerned on the Fourier analysis of the emitted radiation (for 
a review see Ruffini ||) and not on the motion of the particle. Along the lines 
of radiation emission analysis, more recent work has been mainly developed in 
Japan (Osaka and Tokio), United States (Salt Lake City and Caltech), while 
symbolic computing on perturbations has been pursued in Canada (Kingston) 
and Italy (Roma, La Sapienza, ICRA). 

Let us examine the known approaches to solutions of the RW equations: a 
Fourier antitrasform from an approximate frequency solution (where low or high 
frequencies are neglected, and thus also the corresponding antitransformed and 
yet unspecified time-parts of the trajectory of the incoming test mass) presents 
controversial issues of interpretation; a numerical time solution is not adequate 
for comprehension of the different contributions in the geodesic equations. 
Therefore, an analytic solution, in time and obviously radial coordinate do- 
mains, remains the only way for studying the problem of motion. 

Past work was confined to a particle falling in an unperturbed Schwarzschild 
metric, but radiation reaction requires a geodesic path through the emission of 
radiation, in a perturbed Schwarzschild metric. 

This novel approach for calculation of the motion and radiation reaction for 
the two-body problem (body plus particle, the small parameter m/M being 
the ratio of the masses) is based on the following concept: in the background 
curvature given by the Schwarzschild geometry rippled by gravitational waves, 
the geodesic equations insure the presence of radiation reaction also for high 
velocities and strong field (the expansion parameter being m/M and not v/c 
or 4>/c). This concept is applicable to any orbit, but radial fall is of interest 
due to the non-adiabatic regime (equality of radiation reaction and fall time 
scales), in which the particle has to locally and immediately react to the emitted 
radiation. The energy balance hypothesis (emitted radiation equal to the vari- 
ation in the kinetic energy) is only used for the determination of the 4-velocity 
via the Lagrangian and for the normalization of divergencies; the solution in 
time domain of the Regge-Wheeler-Zerilli-Moncrief radial wave equation deter- 
mines the metric tensor expressing the polar perturbations, in terms of which 
the geodesic equations are written; finally radiation reaction shall be identified 
by subtraction of the Oth-order terms in the geodesic equations [^, Q . 
The long-term aim of this work is the identification of radiation reaction in the 
trajectory up to the horizon, without the assumption of adiabacity and with 
minimal use of the energy balance postulate. 

Future developments may include: i) a general method for the determination of 
the motion of small objects in any orbit; ii) a contribution to the solution of the 
inhomogeneous 2nd-order equations, which energy- momentum tensor is based 
on the geodesic equations with perturbation effects and yet to be found; iii) a 
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post-Schwarzschildian formalism. 

In this paper we are concerned with the partial fulfillment of the first step, 
namely the solution of the associated homogeneous equation. 

Radiation reaction is a fundamental concept in bodies motion theory, but 
also has relevant implications on detector's templates since the capture of stars 
by black holes is a source of gravitational waves. 

For a brief critical review of the status of comprehension of radiation reaction 
and references, please refer to HQ]. 

2 The equation 

The RWZM equation for polar perturbations is: 



dr* 2 dt 2 

where 



-V l (r)9i(r,t) = St(r,t) (1) 



r = r 



2Mln (^T7 _1 ) ( 2 ) 



„2M 

is the tortoise coordinate and the potential Vi (r) is 



. . / 2M\ 2A 2 (A + l)r 3 + 6A 2 Mr 2 + 18AA/ 2 r + 18M 3 
W=(, 1 -— J r 3 (Ar + 3M) 2 (3) 

Further A = \{l — + 2) and Si(r, t) is the 2'-pole source component and, 
for a radially falling particle, is expressed by: 



M 6Mr 
r \r + 3M 



where k = 4m y/ (21 + l)ir and ro(t) (known geodesic equation of motion in 
unperturbed Schwarzs child metric) is the inverse of: 



S[r - r (t)] 
(4) 



/ r \ 12 AM ( r \ 3 2 
\2MJ 3 \2M J 



-i 



(5) 

Eq. U reveals that the time domain spans from — oo to +oo. The total energy 
radiated after a given to per / mode is: 



3 



C*>"* (6) 



64tt (I - 2) 

3 The analytic solution 

Eq. (|l]) can be rewritten in terms of t and r* (tortoise coordinate) only, resulting 
into a partial differential equation with constant coefficients but solely via an 
approximate inverse function r(r*) of eq. (0) and thus resulting into an ap- 
proximate p.d.e. Further, the solution is most interesting at r* — —oo where 
a detailed analytic study is not possible. In the (t,r) domain instead, eq.(0) 
becomes: 



1 d 2 ^ t 1 - A{r) (PS>i d 2 ^i 



A 2 {r) dr 2 rA 2 {r) dr dt 2 



Vi(r)*i = Si(r,t) (7) 



dr* r 

where A(r) = —— = . The initial value problem is well defined (ty 

w dr r-2M y 

and '5 are zero, i.e. particle at rest at infinity as the boundary conditions which 

affirm that the same functions are zero at infinity at all times). The associated 

homogeneous equation of eq. (H) can be reduced to an ordinary differential 

equation with the position \l/(r, t) = R(r)T(t). Dividing eq. @) for RT it is 

obtained that (dropping the 1 notation): 

1 d 2 R 1 - A dR 1 d 2 T 

R~A 2 ~d~r 2 ~rRA 2 ~d~r ~ T~dt 2 ^ ~ ^ ( 8 ) 

The first, second and fourth term are not dependent upon t and form time 
constants. Thus also the third term must be a constant in time and eq. (||) is 
autonomous M: 

-— =r> (9) 

T dt 2 [ ' 

Let us pose dT /dt = u and hence d 2 T/dt = uut where the lower subscript 
indicates a derivative to T. Eq. (EJ) becomes: 

^UU T = T 2 (fO) 

which solution is T = e ±r . Eq. (||) becomes: 
^d*R_l-AdR_ 

RA 2 dr 2 rRA 2 dr 1 ' 

or showing the dependence of A(r) with r: 

d 2 R 2M dR r 2 R , , 

;(r 2 +V) = (12) 



dr 2 r{r - 2M) dr (r - 2M) 2 
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For searching the solution, it is better to recast eq. ( |l2] ) as 
d 2 R 2M dR 



^-^2M)^-^ )R = Q (13) 
where f(r) is given by: 

AVV 6 4^Ar 2 Mr 5 +(9r 2 M 2 +2A 3 +2A 2 )r 4 +2A 2 M(l-2A)r 3 46AM 2 (3-2A)r 2 +18M 3 (l-2A)r-36M 4 

r 2 (r -2M) 2 (Xr + 3M) 2 

(14) 



With the position p = r — 2M, eq. (|13|) is finally recast as: 

{p+ 2rn?[X {p+ 2M)+m 2d ^+ 2M{P + 2mm ^ = 

dp z p dp p A 

(15) 

where 

j(p) = - [X 2 T 2 (p + 2Mf + 6Xr 2 M(p + 2Mf + (9r 2 A/ 2 + 2A 3 + 2A 2 )(p + 2M) 4 

+2A 2 M(l-2A)(p+2M) 3 +6AA/ 2 (3-2A)( / 9+2Af) 2 +18Af 3 (l-2A)(p+2M)-36M 4 ] 

_ _ (16) 

Eq. (|15| ) may be solved with the Frobenius method M. Let us rewrite eq. 
© as: 



a (p)—n^ + — — i^ + -^ R = ( 17 ) 



d 2 R{p) | f3{p)dR{ P ) | 7(g) 
c£p 2 p dp p 2 

where the expressions of a(p) and /3(p) are straightforwardly derived from 
15"|). The functions a, /3 and 7 are regular about the expansion point p = and 
have the form: 



n— n=0 n— 

The indicial equation (aos 2 + ((3q — ao)s + 70 — 0) is: 



4M 4 (2A + 3) 2 s 2 + [4M 4 (2A + 3) 2 -4M 4 (2A + 3) 2 )s-16 J ftT 2 M 6 (2A + 3) 2 = (19) 

and reveals that the roots are S1.2 = ±2KM and differ by an integer (4KM). 
One solution R a (p) is a series of the form: 

00 

R a ( P )=P s J2 R °P n ( 2 °) 

n=0 
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while the other is of the form: 



Rb(p) = lnp 



J2(s + 2KM)R n (s)p n 



- b=-1KM n=0 

(21) 



i?o has an arbitrary value and all other terms depend upon it, e.g.: 

A M 3 [(-16A 2 - 102A - 36)K 2 M 2 + (8A 2 + 12A + 9)KM + (4A 2 + 10A + 9)] 

Rl - ~ 4 akmTi 



(22) 
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